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ABSTRACT 

Spectra of Seyfert Is are commonly modelled as emission from an X-ray illuminated 
flat accretion disc orbiting a central black hole. This provides both a reprocessed and 
direct component of the X-ray emission as required by observations of individual objects 
and possibly a fraction of the cosmological X-ray background. There is some observational 
motivation to at least consider the role that an effectively concave disc surface might play: 

(1) a reprocessed fraction > 1/2 in some Seyferts and possibly in the X-ray background, and 

(2) the commonality of a sharp iron line peak for Seyferts at 6.4KeV despite a dependence 
of peak location on inclination angle for flat disc models. Here it is shown that a concave 
disc may not only provide a larger total fraction of reprocessed photons, but can also 
reprocess a much larger fraction of photons in its outer regions when compared to a flat 
disc. This reduces the sensitivity of the 6.4KeV peak location to the inner disc inclination 
angle because the outer regions are less affected by Doppler and gravitational effects. If the 
X-ray source is isotropic, the reprocessed fraction is directly determined by the concavity. 
If the X-ray source is anisotropic, the location of iron line peak can still be determined by 
concavity but the total reflected fraction need not be as large as for the isotropic emitter 
case. The geometric calculations herein are applicable to general accretion disc systems 
illuminated from the center. 

Key Words: accretion, accretion discs; line: profiles; galaxies: active; X-rays: galax- 
ies; Xrays: stars 
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1. Introduction 



Accretion discs are the standard paradigm to explain a wide variety of luminous galac- 
tic and extra-galactic accreting sources. Steady accretion requires dissipation of gravita- 
tional energy which produces the observed luminosity. Emission from X-ray binaries and 
active galactic nuclei (AGN) is thought to result from accretion onto a central massive 
black hole (e.g. Pringle 1981; Rees, 1984). The X-rays originate from the inner-most 
regions of the accretion flow and probe the associated dynamics and geometry. 

X-ray spectra of Seyfert I AGN have been modeled by a combination of direct and 
reprocessed emission (see Mushotzsky et al. (1993) for a review). The direct component is 
from hot 10 9 K electrons of optical depth < 1. The reprocessing component (e.g. Guilbert 
& Rees 1988; George & Fabian 1991) is composed of: 1) a Compton thick material in 
a moderate state of ionisation believed to be the thin accretion disc at < 10 5 K, which 
produces iron fluorescence features and 2) possibly a Compton thin, highly ionised "warm 
absorber" which produces absorption features below 3 keV. If the disc extends to the inner 
stable orbit, a broad gravitationally red-shifted ~ 6.4 KeV iron K a fluorescence line can 
be produced. Material further inside can also contribute to the line (Reynolds & Begelman 
1997) making a distinction between Kerr or Schwarzchild holes difficult (though see Young 
et al. (1998)). Generally, the iron line shape provides a diagnostic for strong gravity 
(Fabian et al., 1995; Tanaka et al 1996). Similar lines have also been seen in galactic black 
hole candidates. (Fabian et al. 1989; Done et al. 1992; Ebisawa 1996). 

The best studied (post ASCA) iron line is that of MCG-6-30-15 (Tanaka et al 1995; 
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Iwasawa et al. 1996). Its profile varies with the continuum, but has been successfully 
modeled by reflection off of a flat Keplerian disc inclined ~ 30 degrees to the line of sight. 
However, ASCA has observed the ~ 6.4 KeV Iron lines in 22 Seyfert Is (Nandra et al. 
1997) and as a population, they have a dispersion of only ±3 degrees around 29 degrees in 
the predicted inclination angle when modelled with a flat disc. The peak of the line always 
appears at or near the rest frame 6.4KeV. Because the peak location is sensitive to the 
disc inclination angle in flat disc models (Laor 1991), the low dispersion is not expected 
for a distribution of randomly oriented discs — even in the presence of a dusty torus at 
large radii which narrows the angular range over which Seyfert Is, by class, are selected. 
The narrow iron line core near 6.4 KeV correlates with the intensity of the continuum flux 
when averaged over ~ few x 10 4 seconds (Iwasawa et al. 1996). Perhaps this indicates a 
corresponding distance between the reprocessing region and the direct X-ray source. 

A second issue is the ratio of total reprocessed to direct emission. In some Seyferts like 
MCG-6-30-15 this may slightly exceed 1 (Lee et al. 1998; Guainazzi et al. 1999), although 
for a population of 11 objects, the ratio seems to hover around 1 (Matt 1998). Some models 
of the cosmological X-ray background (c.f. Fabian 1992) also suggest that the reprocessed 
X-ray component may exceed the direct component by a ratio > 5. If true, this might be 
explained by flat disc geometries by employing an anisotropic direct X-ray source through 
the inverse Compton process, (Ghisellini et al 1990; Rogers 1991), direct acceleration of 
electrons toward the disc (Field & Rogers 1993), source motion (Reynolds & Fabian 1997), 
or general relativistic (GR) effects (Martocchia & Matt 1996). But other possibilities, such 
as concavity, deserve investigation. It important to note however, that recent models of the 
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X-ray background (Comastri et al. 1995) do not require the reprocessed to direct emission 
ratio much different from 1. 

A disc whose outer parts are thicker or concave with increasing height at larger radii 
may play a role in explaining the ubiquity of the iron line peak at 6.4KeV even if the total 
reprocessed to direct emission ratio inferred for a given object is of order 1. The concavity 
allows a larger fraction of the reprocessed emission to be reprocessed in the outer parts 
of the disc, away from the influence of the doppler and gravitational effects. The total 
reprocessed to direct fraction depends not only on the concavity but also on how isotropic 
the X-ray source is. For an isotropic X-ray source, the total reprocessed fraction will be 
large if concavity plays a strong role in determing where the iron line peak is. For an 
anistropic souce however, the concavity may play a strong role in determining where the 
iron line source is even if the total reprocessed fraction is modest. 

Additional "concavity" beyond that of the simplest Shakura-Sunayev (Shakura 1973) 
discs can reflect a thickening disc with a flare due to the particular vertical structure and 
temperature profile (e.g. Keynon & Hartmann 1987 in the context of stellar discs.) Also, 
azimuthally dependent concavity can result from warping which may be tidally or radia- 
tively driven (Terquem & Bernout 1993,1996; Pringle 1996,1997), or induced by a wind. 
Alternatively, the discs may incur a transition from a thick torus to a thin disc inward. 
Here I do not consider the dynamics in detail and just parameterize the concavity to high- 
light some simple effects related to the above observations. Concavity was considered by 
Matt et al. (1991) primarily for its effect on the iron line equivalent width, but this may 
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not be its most important role. In section 2, I explicitly derive the ratio of reprocessed 
to direct emission as a function of the curvature. Splitting the reprocessed contribution 
into components emanating from inside and outside a critical radius r c , I then derive the 
ratio of contributions to the iron line from the outer (taken to be the "narrow component" 
around 6.4KeV) and inner regions (taken to be the "broad component"). The calculation 
results and line profile examples are discussed in section 3 and section 4 is the conclusion. 
2. Concavity and the Reprocessed Emission 

2.1 Basic Considerations 

Here I take he direct X-ray source is taken to be an isotropic point emitter located 
above a Schwarzchild hole. (Later I will comment on the possibility of an anisotropic 
source.) The exact location of the X-ray source(s) is unknown. To marginally avoid 
considering reflection from material free-falling inside of radii r = 6 (Reynolds & Begelman 
1997), I take the source height to be H e = 10, where the r and H e are in gravitational 
units of R g = GM/c 2 , and M is the hole mass. Incident hard X-rays impinge onto the 
accretion disc and are scattered off an optically thin outer layer (see Matt et al. 1996). The 
X-ray photons also excite fluorescent photons from iron atoms. Below, the 6.4 KeV iron 
line flux is taken to be proportional to the incident flux. This assumes that the ionisation 
parameter is low enough over all r > 6 to ensure cold iron line emission (Matt et al. 1996). 
The calculations then reduce to the geometry of Fig 1. 

2.2 Total Reprocessed vs. Direct Flux Ratios for Concave discs 
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First consider the total reprocessed flux. The ratio of the observed flux in the repro- 
cessed component to that in the direct component is the ratio of the number of photons 
which intercept the disc over the number of which escape directly. For the parameters 
chosen above, the GR effects are sub-dominant in all regions of the disc and comprise a 
maximum 30% correction. Note however, that GR leads to more flux impinging on the 
disc, so ignoring GR at first highlights a different way to achieve a a high reprocessed 
fraction, and thus an overall lower limit. 

In the Euclidean regime, the flux ratio of reprocessed to direct emission is given by 
the ratio of solid angle subtended by the disc to that subtended by the free space to the 
observer. The azimuthal angle drops out by symmetry. Take the disc's reprocessing layer 
to have a height 

h(r) = a(r/6) b [for a(r/6) b < r] ; r [for a(r/6) 6 > r], (I) 

where the curvature index b implies concavity for b > I, and a = h(ri n )/ri n at the inner 
disc radius r in = 6. The two regimes of (I) are imposed to enforce h < r at all radii. The 
dashed curves of Fig. 2 show the b for which h(r out ) = r for different outer radii r out . The 
ratio of reprocessed to direct flux, using the angles shown in Fig. I, is then 

rTv-Tan- 1 [r in /H e ] SiflOdO 

rp I Tp <Jn/2—Tan- 1 [(H e — h(r out ))/r ou t] /r>\ 

frep/^dir ~ ^/2-Tan-^[(H e -h(r mt ))/r mt ] g^Q^ 1 ' 

For a strictly flat disc, the ratio is given by (2) with h{r) = 0, and then F rep /Fdi r < 1- 
2.3 Concavity and the Reprocessed Iron Line 

Consider two components to the reprocessed iron line, motivated partly by the ap- 
proach employed for the best-studied MCG-6-30-15 example (e.g. Iwasawa 1996). Take 
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the first component to peak at the rest frame frequency and the second to be the remaining 
broad line. After choosing a core width, one can derive a corresponding critical radius, r c , 
outside of which all the narrow core emission emanates. The role of concavity becomes 
apparent by comparing the flux ratios from the two disc regions for flat vs. curved discs. 

Iwasawa et al (1996) and Nandra et al. (1997) consider a narrow component width 
~ ±2.5% around the 6.4 KeV peak rest energy in MCG-6-30-15 and for the population of 
22 objects respectively. We can estimate r c by assuming that the spread corresponds to a 
Doppler width enveloping the largest and smallest frequencies at r c . The Doppler shift at 
large r is 

i//i/ e ~(l±Av/c), (3) 

where v is the frequency, v e is the rest frequency = 6.4keV, and Av is the maximum spread 
in velocity. For Av we take the Keplerian speed ~ c^l/r) 1 / 2 , Thus 

r c ~(Az//z/ e )" 2 , (4) 

where Av is the frequency half-width of the core. For Av = 0.025, r c = 1600, while for 
Av = 0.05, r c = 400. 

Using Fig. 1 to calculate the ratio of flux emanating from outside r c to that from 
within r c we obtain 

r-Tr/2-Tan- 1 [(H e -h(r c ))/r c ] ginOdd 

rp / rp _ J-n-/2-Tan- 1 {(H e -h(r ou t))/r out ] /r\ 

-Tout/ -Tin — r - K -Tan- 1 [r in /H s ] . ' ^ ' 

L/2-Tan-^\{H £ -h{r c ))/r c \ Smddd 

Notice from Fig. 1 that the angle bounds in the integral are the same regardless of whether 
h(r c ) > H e or h(r c ) < H e . For a flat disc, (5) is found from taking h = for all r. 
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2.3 Local Approach 

The above approach is the simplest for the Euclidean case, but to include GR and 
actually compute line profiles, GR corrections are more easily incorporated in a formalism 
which integrates over r. An equivalent way to compute (5) is to note that the reprocessed 
line flux from an element of area dA is proportional to the impinging flux f(r,H e ,h) at 
that radius projected onto the area dA (see Fig. 1) 

dF oc g r (r)f(r, H e , h)Cos\ r(dr 2 + dh 2 ) 1/2 = g r (r)f(r, H e , h)Cos\{\ + dh 2 /dr 2 ) 1/2 rdr 

oc g r (r)(r 2 + (H e - h) 2 )' 1 ^ + dh 2 j dr 2 ) x l 2 rdr, (6) 

where f(r,H e , h) cx (r 2 + (H e — /i) 2 ) -1 . The angle A is between the normal to the area 
element and the X-ray source direction (Fig. 1) so that 

CosX = Cos(n/2 - Tan'^dh/dr) + Tan- 1 [(/i(r) - H e )/r}). (7) 

The quantity g r (r) = (1 — 2/r)~ 4 (l — 2/r) 1 / 2 (l + z)~^ is the product of three correction 
factors to an otherwise Euclidean formula: The first is the Doppler + GR correction to the 
Euclidean illumination function as approximated by Reynolds & Begelman (1996). The 
second is the correction to the area measure in the integrand. The third factor (1 + z)~ 4 
is the red-shift correction. This is given by (1 + z) = {j>^u^) ern / E b s where subscript em 
stands for emitted location (i.e. disc), E ohs is the photon energy measured by the observer, 
and p^Ufj, is the product of the photon 4-momentum and bulk motion 4- velocity along the 
sight line. At finite inclination angle, some of the photons from the inner regions populate 
the narrow core. In addition, some inner region photons may be lost as they impinge on 
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the outer disc and incur a second reprocessing. The formulae for a face on inner disc thus 
provide a lower limit to the ratio F out /F in for all inclination angles. From Reynolds & 
Begelman (1996), we then obtain (1 + z)~ 4 = (1 — 3/r) 2 . Notice that the illumination 
function correction and the red-shift correction compete. The former enhances, while the 
latter decreases the flux. The net effect is an increase in the redshifted component and an 
increase flux from the inner regions, and thus a decrease in (5). I ignore the tiny corrections 
to g r (r) which result from a finite h. 

Using (7) in (6) gives 

[XT- g r (r')dr'r'(l + ( fr) 2 ) l ' 2 (r' 2 + (H e - hf)^Cos\(r')] 
° Ut/ ^ Uri 9r(r')dr'r'(l + (f^W 2 + (#e - h)*)-iCos\(r')] ' 

For a flat disc, dh/dr = and then CosX = (H e - h)/(r 2 + (H e -h) 2 ) 1/2 . In the limit 
g(r) = 1 for all r, the results from using (8) are identical to those from (5). In the next 
section I discuss the results of the above ratios and show some line profiles. 

3. Results, Line Profiles, and Discussion 

The solid curves of Fig. 2 show F rep /F dir as a function of b for three values of r out: 
using a = 1/50, H e = 10, and b > 1. The restriction to concave curvature (i.e. b > 1) 
means that the entire disc surface sees the X-rays. The imposed restriction that h < r 
leads to a maximum ratio of 5.4 in Fig 2. For the above parameters, a Euclidean flat disc 
would give only F rep /Fdi r ~ 0.87 for an isotropic source. Note that if the X-ray source 
were suitably anisotropic, the ratio F rep / ' F^ir need not be significantly greater than 1 even 
for the optimal b. In that case, the Euclidean flat disc ratio for the same anisotropy would 
still be of order 5 times less for the optimal b. 



Fig. 3 shows F out /F in , for different values of r out and r c . (To be conservative, Fig. 3 
employs (8) in order to include the rough GR corrections which lower the curves relative 
to the Euclidean case, making the effect ~ 30% less pronounced. By contrast, for Fig 2, 
the conservative lower limit is given by the Euclidean equation (2)). There are two regions 
for each curve of Fig. 3. For low 6, the gain from orientation toward the source wins 
over the decrease in flux from the extra distance to the disc at a given r and the curves 
rise. However since H e « r c , above the b at which h(r)/r = 1, the increasing distance 
between the X-ray point source and the disc for larger r wins and the curves then decline. 
A Shakura-Sunaeyev (Shakura & Sunaeyev 1973) disc corresponds to b = 9/8. For a flat 
disc the ratios are low: for r c = 400, F out /Fi n ~ 0.03 and for r c = 1600, F out /Fi n ~ 0.007. 

The r c = 400 solid curve lies completely above the r c = 1600 solid curve in Fig. 3 since 
r out is the same in both cases. This contrasts the dashed curves which address a different 
issue. The approximate 2 — 4 x 10 4 sec delay between an increase in continuum emission 
and the response of the narrow core line for e.g. MCG-6-30-15, (Iwasawa et al. 1996) 
motivates considering that some reprocessing material resides at the associated distance 
from the X-ray source. For a 1O 7 M hole, this corresponds to a distance of 400 — 800R g . 
The Keplerian velocities at this r c are consistent with a few percent frequency width of 
the core. This motivates determining the amount of reprocessed emission coming from 
r c < r < 10r c and the results are the dashed curves of Fig 3. Because of the narrow range 
in r, the curves are lower than for the solid curves but there is still a range of b for which 
such an outer "ring" can produce F out /F in > 1/3, even for r out = 4000, which is of order 
that required for MCG-6-30-15 (Iwasawa et al. 1996). By contrast, for flat discs from (2), 
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even the more favorable case (r out = 10r c , r c = 400) gives only F out /F in = 0.03. The 
dashed curves in Fig. 3 cross because a decade in r c for higher r c is larger, but farther 
out, than a decade in r c for lower r c . 

The fact that emission in the narrow peak could originate from large radii means 
that the ubiquity of observed peaks near 6.4keV in the 22 Seyferts studied by Nandra et al 
(1997) would not be as sensitive to the inner disc inclination angle. The maximum velocity 
dispersions at such large radii are small (e.g. Au < 5% at r = 400). For a strictly flat 
disc, some tuning in inclination angle is required to produce the location of the observed 
peak because gravitational+transverse Doppler red, and blue shifts conspire to produce 
the particular peak values. Figure 4 shows some line profiles using the formulae of Fabian 
et al. (1989) (Esin 1998) but with a modified emissivity function to include the flat and 
concave disc cases. The emissivity function of Fabian et al. (1989) was e oc (r in /r) 2 and I 
replace this by the factor in (8): (1 + (^) 2 ) 1/2 (r 2 + (H e - h) 2 )- l Cos\{r). This simplifies 
for a flat disc as discussed below (8). As expected from Figs. 2 and 3, Fig 4 shows the 
narrow peak at 6.4KeV for the value of b = 1.7 even at an inclination angle of 40deg. The 
effect weakens significantly for b = 1.3 as expected. Secondary reprocessing (Matt et al. 
1991) of inner disc radiation by outer disc radiation is not considered. This would be less 
likely to effect the peak because the peak is produced from emission at large radii. Also, 
the constraint h < r means angles of inclination < 45deg are less affected by shadowing. 

The above results show that for a range of b and r out: the ratios F rep /F dir and F out /F in 
can be more than an order of magnitude larger than for flat discs. The concavity may 
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thus predominantly affect the total reprocessed fraction or the line shape rather than the 
line equivalent width (Matt et al. 1991). The results are not strongly sensitive to the 
parameters a, H e , or ri n for H e ,r e > 6, but for ri n ,H e < 6 a more detailed inclusion of 
the ionisation fraction and shadowing is required. 

The value of b which provides both the largest F out /F in and F rep /Fdi r for the range of 
r out considered is 1.5 < b < 1.9 out to the radius for which h(r)/r = 1. (The value b ~ 1.5 
is that of an isothermal disc.) A disc which changes from a thin to thick disc/torus at 
r > r c and has reprocessng material in the torus, may be approximated by (1). Some excess 
in reprocessed emission can also result from a warped disc, possibly tidally or radiatively 
driven (e.g. Terquem & Bertout 1993,1996; Pringle 1996,1997) or induced by wind torques, 
but the azimuthal dependence must then be considered (e.g. Terquem & Bertout 1993). 
The role of an anisotropic X-ray source would change the total reprocessed fraction, but 
not the influence on the line peak location, or the comparisons to a flat disc. 

4. Conclusions 

Some simple but pronounced effects of a concavely curved accretion disc are captured 
by estimating the total reprocessed vs. direct flux and the relative flux emanating from 
inside and outside of a critical radius r c . For a range of 6, the reprocessed flux from r > 400 
and even r > 1600 can contribute significantly to a rest frame iron line peak even when 
the (X-ray) radiation source is at located at H e ~ 10. This alleviates some sensitivity 
of disc model predictions (e.g. Laor 1991) to the disc inclination angle, as seen in Fig 4. 
Reprocessing at large distances also predicts a time delay between changes in the direct 



12 



continuum emission and the reprocessed emission which may be observed in some Seyferts 
(Iwasawa et al. 1996). 

In addition, concave discs may produce F rep /Fdi r > 5 for an isotropic X-ray source 
which provides another means for some AGN models of the X-ray background to account 
for the high required reprocessed fraction (c.f. Fabian 1992). However, if the X-ray back- 
round and Seyferts do not generically require a high reprocessed fraction (Comastri et al 
1995; Matt 1998), the ubiquity of the iron line peak at 6.4 KeV could still be influenced 
by concavity when compared to a flat disc: The concavity affects not only the total repro- 
cessed fraction for a given X-ray source anisotropy, but the relative reprocessed fraction 
from different parts of the disc. Finally, note that the disc could be curved or could change 
from a thin to a thick disc with entrained reprocessing material, such that the effective 
global reprocessing geometry is approximated by a concavity. 
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Figure 1 Caption: Schematic of X-ray source (small circle) above disc surface (thick 
black curve). Two values of 6(r) and A(r) angles are shown for illustration as used in (2), 
(5) and (8). Regime (a) H e < h, and (b) H e > h. 

Figure 2 Caption: Plot of h(r)/r (dashed curves) and F rep /F dir (solid curves) vs. 
b for r ou t = 10 5 , 16000, 4000 from left to right in each curve group. Values of H e = 10, 
ri n = 6, and a = 1/50 are used in all curves. The flattening is the result of imposing 
h(r) < r for all r in (1). The X-ray source is assumed to be isotropic for this graph. 
Simple anisotropies could be included by multiplying the y-axis by a constant fraction. 

Figure 3 Caption: F out /F in vs. b from (8) for r out = 10 5 (solid curves) and r out = 
10r c (dashed curves). The top solid and dashed curves as measured at b = 1 have r c = 400, 
and the bottom curves have r c = 1600. The dashed curves cross and merge with the solid 
curves above the b for which h(r)/r = 1. The condition h(r) < r is imposed, but the 
presence of the down-turns only requires r c » Ti n . 

Figure 4 Caption: Line profiles for flat and concave models based on Fabian et al. 
(1989) with r out = 10 4 with emissivity function modified as discussed in text. Solid lines 
are for flat discs at inclination angles of 40 deg (broader curve) and 30 deg (narrower curve) 
respectively. The concave disc curves match onto the wings of the flat curves. There are 
4 curved disc profiles signatured by the height of their peaks. The largest has b = 1.7 and 
30 deg inclination. The next largest is b = 1.7 and 40 deg. The next is b = 1.3 and 30 deg, 
and the lowest is b = 1.3 and 40 deg. 
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